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Abstract 



We have carried out extensive simulations of the angle-resolved photoe- 
mission (ARPES) intensity in Bi2212 within the one-step and three-step type 
models using a first-principles band theory framework. The focus is on under- 
standing the behavior of emissions from the antibonding and bonding bands 
arising from the Cu02 bilayers around the M{7r, 0) symmetry point. The 
specific issues addressed include: Dependencies of the photointensity on the 
energy and polarization of the incident light; character of the initial and fi- 
nal states involved as well as the spectrum of the relevant final states; and, 
changes in the spectral intensity as a function of the perpendicular compo- 
nent, k±, of the momentum of the photoelectron. Considerable insight into the 
nature of individual transitions is adduced by examining the momentum ma- 
trix element for bulk transitions within the solid and by further decomposing 
this matrix element into contributions from various atomic sites and angular 
momentum channels. These results indicate that, via remarkable interference 
effects, the ARPES matrix element can in particular cases help zoom in on the 
properties of electrons excited from specific sites and/or angular momentum 
channels even in a complex material. 
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I. INTRODUCTION 



It is well recognized that the photoemission process involves an interplay between the 
bulk and surface phenomena as the excited electron must be ejected from the surface in order 
to reach the detector. The intensity measured in an ARPES (angle-resolved photoemission 
spectroscopy) experiment, therefore, fundamentally involves the properties of the relevant 
initial and final states in the presence of the surface. As a result, the ARPES intensity 
can differ greatly from the spectral density of the initial state involved in the excitation. 
An understanding of this difference — the so-called 'matrix element effect', is essential for 
a satisfactory interpretation of the ARPES spectra, and this is particularly so in complex 
materials such as the high-Tc's 

Bearing these considerations in mind, we have carried out extensive simulations of the 
ARPES spectra in Bi2Sr2CaCu20s (Bi2212, BISCO) within the first-principles band the- 
ory framework. The pristine phase with the body centered tetragonal lattice structure is 
assumed; the computationally more demanding case of the orthorhombic and modulated 
lattices is hoped to be taken up in the future. We focus on the excitation of states in the 
vicinity of the Fermi energy (Ep) around the M (vr, 0) symmetry point which have been the 
subject of intense scrutiny and controversy in the high-Tc literature. For example, the issue 
of bilayer splitting— expected to be relatively large within the conventional LDA-based band 
theory picture at M in BISCO, has been debated for some time, even though a consensus in 
favor of the existence of a bilayer splitting appears to be emerging in the recent months 0-§| . 
The fundamental question has been whether or not strong correlations localize electrons in 
the Cu02 layers, yielding a 2D electron gas and a zero bilayer splitting with significant 
implications for the mechanism of superconductivity in the cuprates. This study gives 
insight into these and related aspects of the ground state electronic structure by delineating 
the connection between the ARPES spectra and the character of the underlying initial and 
final states. 

For this purpose, we have developed the one-step methodology for modeling the 
photoemission process so that arbitrarily complex systems with many basis atoms can be 
treated P|,^,|9|-pn| . The previous one-step work has been limited to a maximum of two atoms 
per unit cell [0 which is not adequate for the cuprates; BISCO, for example, involves 30 
atoms per conventional unit cell even for the pristine lattice. In the one-step approach, 
the interaction of light with the solid is treated as a single quantum mechanical event and 
the artifical distinctions between the processes of excitation, transport and ejection of the 
photoelectron invoked in the earlier three-step models p!3HT7| are not made. The one-step 
scheme is thus inherently more satisfactory than the three-step model since the surface is 
incorporated in the latter in an ad hoc manner. On the other hand, it often tends to be 
difficult to identify the behavior of individual transitions within the one-step methodology 
where the initial and final states must be damped, causing these states to broaden and 
overlap, blurring distinctions between various transitions. For this reason, we have also 
developed the three-step model for complex lattices [T^. Useful insight is adduced here into 
different contributions to the ARPES spectrum of BISCO by examining the momentum 
matrix element for bulk transitions in the solid in the spirit of the three-step approach. 

The specifics of some of the main issues addressed in this article and highlights of the new 
results presented are as follows. The dependencies of the ARPES intensity on the energy and 
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polarization of the incident light are clarified for the excitation of states near the M point. 
The response of the antibonding and bonding combinations of states from the Cu02 bilayers 
is found to display a strong energy and polarization dependence. Moreover, the behavior 
of the antibonding and bonding states differs significantly so that by a judicious choice of 
these experimental parameters even closely placed bilayer bands could be distinguished. 
Also, photon energies where the intrinsic cross-section for exciting certain states of the 
pristine lattice is large will be best suited for studying these states; in contrast, when this 
cross-section is relatively small, the effects of deviations from the tetragonal symmetry of 
the system via various distortions and modulations will become more prominent. Further, 
we show how some systematics of the polarization dependencies can be understood in the 
cuprates via straightforward arguments involving the symmetries of the initial and final 
states. 

Turning to questions related to individual transitions, we discuss the nature of the mo- 
mentum matrix element for bulk excitations within a three-step type model in terms of 
contributions from various basis atoms and different angular momentum channels. The 
character of the bonding and antibonding initial states as well as that of the relevant final 
states at M is clarified. An analysis of these results reveals dramatically how different con- 
tributions to the momentum matrix element can interfere constructively or destructively to 
enhance or supress the weights of photoelectrons excited from particular atomic sites and/or 
specific angular momentum channels. This remarkable result which, to our knowledge, has 
not been recognized previously in the literature, hints that the ARPES matrix element may 
allow us to zoom in on the properties of very specific transitions even in complex unit cell 
materials. 

Concerning relevant literature, Refs. |[l|,H,P, |l0l|l8i[l9l discuss some of our earlier work 
on BISCO, Nd2-xCe^CuOi (NCCO), YBa2Cu^0j (YBCO) and Y Ba2CuiOsiYl24:). In 
particular, Ref. m emphasizes the importance of matrix element effects in ARPES spectra, 
while Ref. undertakes an extensive comparison between the one-step predictions over 
a wide region in the momentum space when the initial state energy is held fixed at the 
Ep with the corresponding experimental ARPES results giving the Fermi surface map in 
BISCO. There is of course a large body of available work on the photoemission theory and 
Refs. [pil|-p5| would provide a brief selection. On the experimental side as well the ARPES 
studies of BISCO are far too numerous to be cited with any completeness, and we refer only 
to a few representative papers. |3-|6|,|26|-05[] 

An outline of this article as is follows. The introductory remarks are followed in Section 
II by a consideration of formal matters; a few equations for the one-step model are pre- 
sented (Subsection IIA), and the problem of computing the momentum matrix element in a 
general lattice is formulated (Subsection IIB). Section III delineates the relevant features of 
the lattice geometry (Subsection IIIA) and the band structure of BISCO (Subsection IIIB). 
The results of Sections II and III serve as a source of frequent reference in the subsequent 
discussion. Section IV is divided into several subsections to go over various salient features 
of the ARPES spectra from M. The aspects addressed concern: Photon energy dependence 
(Subsection IVA); Effect of final state width (Subsection IVB); Polarization dependence 
(Subsection IVC); and, the final state band structure and momentum matrix element (Sub- 
section IVD). Section V takes up the question of the nature of the transition matrix element 
and the origin of specific spectral features in terms of the characters of the initial and final 
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states and is divided into three parts dealing with: Contributions of different atomic sites 
(Subsection VA); contributions of various angular momentum channels (Subsection VB); 
and, comments on the k± dependence of the spectral intensities (Subsection VC). Finally, 
Section VI summarizes our results and makes a few concluding remarks. 



II. THEORETICAL CONSIDERATIONS 
A. General aspects 

We begin by recalling that the ARPES intensity in the one-step model can be expressed 



as 



j = -l/m<k|||G+ AG+A^Galkii > (1) 

where the ket |k|| > denotes a free electron state with momentum, ky, parallel to the crystal 
surface, G2 (Gi) is the retarded (+) or advanced (-) one-electron Green function for the final 
(initial) state and 

A = e/i/2mc(p ■ A + A ■ p) (2) 

is the interaction Hamiltonian in terms of the vector potential A of the incident photon field 
and the electron momentum operator p. 

For our purposes, we rewrite Eq. ^ by using the representation 

- -ImGt = y^\i> Bii< i\ (3) 
^ i 

of the Green function Gi as a sum over the spectral functions for various initial states i 
which yields 

I = < kiilG+A^ |z > Bii < tlA^G^m > (4) 

i 

= J2Bii< kylG+Alz >< t\A^ G^lkw > (5) 

i 

= J2B,,\<f\A\z>\' (6) 

i 

Here we have replaced |k|| > by the final state wavefunction |/ >. Substituting form for A 
in Eq. H, and neglecting the term proportional to V- A by assuming that the vector potential 
A varies smoothly, i.e. wavelength of light is much larger than the unit cell dimensions ||36|| , 
we obtain 

I = eh/mcJ2Bii\A- < f\p\i> (7) 

i 

Notably, the summation over i not only encompasses different initial state bands, but also 
involves an integration over k^. A computationally practical formula follows by invoking 
the identity 
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[H, p] = Hp-pH = -ihVV 



(8) 



which is vahd when the potential term is p independent. The matrix element in Eq. |^ is 
easily recast using Eq. || as 

A- < /IpK >= -^^^A- < f\VV\i > (9) 
where Ei and Ef are the initial and final state energies, respectively. 



B. Momentum matrix element 

The nature of the photointensity based on the one-step Eq. |l| can be delineated in terms 
of the momentum matrix element < f\p\i > of Eq. 0. However, the wavefunctions \i > and 
1/ > occuring in Eqs. |^ - ^ are complicated quantities which include multiple scattering ef- 
fects in the presence of the surface and damped initial and final state propagators. Although 
the full one-step calculation should always be kept in mind, insight into the results may be 
obtained by simplifying the situation and replacing \i > and |/ > by the corresponding 
Bloch wavefunctions ip in the bulk crystal and the resulting matrix element < ipflplipi > 
for optical transitions. It should be emphasized that the optical matrix element must in 
general be supplemented with other information specific to the photoemission experiment 
(e.g. symmetry of the final state) before connection with a particular measurement can be 
made. This issue will be addressed as needed in the remainder of this article. 

The formalism for evaluating the optical matrix element for a general lattice within the 



KKR scheme is given in Refs. |1T7| , |37|| ; an outline of aspects relevant for the purposes of this 
article is as follows. The starting point is the KKR wavefunction for a general lattice which 
can be expressed within the unit cell as [0^ , ^] 



m = E ^'C^LR?ir)YL{n) (10) 

where L = (/, m) is a composite angular momentum index and /3 denotes different basis sites. 
Cf are expansion coefficients and -Rf (r) is the radial part of the Bloch wavefunction on site 
j3 in angular momentum channel L. Yl{^1) are real spherical harmonics. Using expressions 
of Eq. |10| for the initial and final states yields 

< i>f\m >= E E ej-'-'c^,:c'Xi'S2,L' (11) 

a,f3 L,L' 

where the primed indices refer to the final state and is a unit vector along the direction 
a. Bill involves an integral over the radial part of the initial and final state wavefunction 
and their derivatives and is given by 
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p^Pi - p^p; 



I' = 1+1 



Tfj 





+ ^ 



(12) 



E,-E 



P'/P - P;P/ 



r=rp 



I' = I 



Here, Pi{r) = rRi{r), rp is the radius of the (5^^ muffin-tin sphere, Ei and Ef have been 
defined in connection with Eq. ^ above, and 







j YL'{VL)Yi^^{VL)YL{VL)dVL V = l±l 



(13) 



Ql j^i are seen to contain angular momentum selection rules for dipole transitions. Yi ^ for 
a values of -1,1 and is proportional to x^y and respectively. 

It is convenient to decompose the sum over j3 on the right hand side of Eq. 11 for the 



^th component as 



M^=<i,f\p^\i,,> =Y.M2 

7 5 

7 5 L,L' 



with 



7,(5 



{cifycrBi^Gl,, 



(14) 
(15) 
(16) 

(17) 



The indices 7 and 6 in Eqs. p!4|-p!q are keyed in with the pristine (tetragonal) crystal structure 
of BISCO detailed in Section IIIA below, and taken together, encompass the sum over all 
atoms (3 in the unit cell in Eq. 10. After the sum over the angular momentum channels L 
and L' in Eq. ^ has been carried out, the index 6 which only takes two values (1 and 2) 
in Eq. |15| sums over the contributions of pairs of atoms (other than the Ca atom) related 
by mirror symmetry of the lattice with respect to the Ca layer. In this sense, 6 may be 
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thought of as a "pairing" index. 7 then takes on eight distinct values which include seven 
pairs of atoms {Bi, Obi, Sr, Osr, Cu, Ocu,x, Ocu,y) and the Ca atom as a "site index". The 
aforementioned indexing is by no means unique and other useful variations can be envisioned 
in the present case and certainly for other lattice types. 

III. RELEVANT FEATURES OF THE LATTICE GEOMETRY AND 
ELECTRONIC STRUCTURE OF BISCO 

A. Lattice structure 

Since various structural features of the lattice are invoked frequently in the discussion, 
it is appropriate that we state our nomenclature clearly. For this purpose, two different 
cross-sections through the conventional body centered tetragonal unit cell, together with 
the detailed arrangement of atoms in different layers, are shown in Fig. 1. The distinct 
atom "pairs" to go with the index 7 in Eqs. are identified in Fig. 1(c). Values of 7 

from 2-8 involve an atom and its mirror partner reflected in the Ca layer, yielding a total of 
15 basis atoms per primitive lattice point or 30 atoms in the conventional unit cell. The Ca 
atom (7 = 1) is seen to lie at the center of symmetry in the upper part of Fig. 1(a). The 
two Cu02 planes sit above and below the Ca layer and consist of the two Cu atoms in the 
basis set denoted by the (7,5) indices (2,1) and (2,2). Ocu,x and Ocu,y are given similarly 
by 7 = 3 and 4, each with 5 = 1,2; note that although Ocu,x and Ocu,y are identical in 
terms of the crystal potential in the tetragonal case, it is still useful to make this distinction 
because states along x and y will in general respond differently depending upon the direction 
of polarization of the incident light. Next come the SrO layers with Sr (7 = 5) and the 
associated apical Osr atoms (7 = 6), followed by the BiO layers involving the Bi (7 = 7) 
and Obi (7 = 8) atoms. 

B. Electronic structure 

Fig. 2 shows the energy bands in BISCO along the symmetry lines F — M — X — F. The 
crystal potential underlying these computations differs somewhat from the selfconsistent 
LDA potential in that the Bi-0 pockets around the M-point have been lifted above the 
Ep in order to account for their absence in the ARPES measurements. The two bands in 
the vicinity of the Ep are related to the CUO2 planes; the splitting between a and h is the 
bilayer splitting arising from interaction between the two CUO2 planes. At the M-point, 
the wavefunction of the upper band a is antibonding or antisymmetric with respect to the 
Ca plane (along the z direction) while that of the lower band b is bonding or symmetric. 
The size of the bilayer splitting is generally larger along F — M compared to the F — X 
direction. We have chosen to work in the present simulations with a potential that yields 
a rather large bilayer splitting because this makes it easier to identify the signatures of the 
bonding and antibonding bands h and a in the theoretical ARPES spectra [0. However, 
computations have been repeated for numerous other cases where the bilayer splitting was 
varied by placing suitable potential barriers between the CUO2 layers, and the results indicate 
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that our conclusions concerning the nature of these states and their behavior in the ARPES 
spectra are robust to such variations. 

The analysis of the recent high resolution ARPES spectra leaves little doubt that the 
Fermi surface of BISCO consists of two distinct sheets One is the standard large 

hole-like sheet centered around the XorY symmetry point, and a second sheet arising from 
a band lying very close the the Ep which could be either slightly hole- or electron-like. The 
precise position of the Ep will depend of course on the doping level of the system. In any 
event, it is often useful to avoid convo luting the computed spectra with the Fermi function 
so that the nature of the spectra over a wider range of ky and energy values, including 
unoccupied states above the Ep, can be elucidated. Notably, the k^ dispersion of bands 
terminating at points a and b, even though it varies with ky, is quite small (order of a few 
meV at M) and is not important in much of our discussion. 



IV. SALIENT FEATURES OF THE ARPES SPECTRA FROM THE M POINT 



A. Photon energy dependence of emission intensity from bonding and antibonding 

states 



Fig. 3 shows the ARPES intensity from the antibonding and bonding states (a and b) 
at the M point as a function of the photon energy. For our illustrative purposes we have 
taken the light to be polarized such that only the x-component of the vector potential A is 



non-zero. The initial state width corresponds to T,'- = 100 meV The final state width is 
set using S" = 1 eV— we return to this point in Section IVB below. The emission intensity is 
seen to change dramatically with photon energy. The bonding state emits strongly around 18 
and 23 eV with a feature around 33 eV, while the antibonding state possesses a somewhat 
broad region of high emissions extending from 14-26 eV. The intensity for both states is 
relatively small from 27-40 eV and below 11 eV. These results have significant implications 
for ARPES experiments. For example, emission maps from the Ep can be expected to look 
quite different at various photon energies; at an energy such as 22 eV, the emissions from 
the pristine lattice will be dominant, while at an energy such as 28 eV, the contributions 
due to modulations of the underlying symmetry (e.g. orthorhombic distortion, superlattice 
modulation, etc.) will generally become more prominent |^. Energies such as 18 or 23 
eV will tend to highlight the bonding state while the antibonding state will likely be more 
intense around 22 eV, effects of final state broadening and other factors notwithstanding. 
We should particularly keep in mind that there are inherent uncertainities of a few eV's in 
locating the final states in the first principles computations, and that the positions of various 
features in Fig. 3 will generally possess similar uncertainty. We have carried out additional 
computations (not shown in the interest of brevity) and find that the aforementioned type of 
variations in intensity which selectively highlight different aspects of the electronic structure 
continue to be manifest at photon energies much higher than the 40 eV upper limit in Fig. 
3. 
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B. Effect of final state width 



Fig. 4 shows how the plot of Fig. 3 for the antibonding initial state changes when the 
final state width varies corresponding to values from 20 meV to 2 eV |^ . The results for 
small Sj make it clear that the spectrum is intrinsically made up of numerous transitions 
with different weights. With increasing these transitions broaden and overlap to produce 
varying spectral shapes of Fig. 4. For H'j- < 100 meV, little change in shape occurs because 
the spectra are controlled by the initial and not the final state width, the former having a 
constant value given by = 100 meV in all cases. Note that to reproduce the observed 
linewidths, used in the computations will need to account not only for the damping 
of the final state reflected via the imaginary part of the self-energy, but also the effects 
of experimental parameters such as the temperature of the measurement and the acceptor 
angle of the analyzer (typically of order of a few tenths of a degree in the high resolution 
ARPES setups). In this sense, is perhaps better viewed as a computational rather than a 
physical parameter in the theory with values ranging from 1-4 eV, where 4 eV is appropriate 
for 40-60 eV photons |]42|j43[| . Incidentally, computations become more demanding as the 
value of decreases because the mean free path of the outgoing electron increases and 
an exponentially larger number of layers must be included in the calculations to obtain a 
converged result. 



C. Polarization dependence 

Insofar as the polarization dependence is concerned, note first that in Fig. 5 light is 
assumed incident normally to the surface and that the detector lies in the x — z plane with 
the X-axis defined by the Cu — bond direction; the polarization vector of light is varied 
holding all other experimental parameters fixed. The inserts in the shape of number eight 
are polar plots of the intensities of various peaks (a, b, c and c') as a function of the azimuthal 
angle between the polarization vector and the x-axis. The intensity is maximum for all 
peaks at = 0, i.e. when the light is polarized along the Cu — bonds. Under a 90° 
rotation the intensity nearly vanishes as the polarization vector turns perpendicular to the 
plane of the detector. 

Insight into the polarization dependence discussed in the preceding paragraph can be 
gained by examining the behavior of Eq. |^ for the geometrical setup of Fig. 5. Since the 
detector lies in the x — z plane of mirror symmetry, the final state \f > must obviously 
possess even symmetry with respect to this plane in order to be observable. On the other 
hand, or (or, equivalently, dV/dx or dV/dy from Eq. ^ is even but Py is odd under 
reflections in this mirror plane. It follows then that only the matrix elements < f\px\i > and 
< /bzK > can be non-zero for even initial states while for odd initial states only < f\py\i > 
can be non-zero. We now express A in spherical polar coordinates as 

A = AqC = y4o(cos (f) COS 6e^ + sin cos Osy + sin Oe^) (18) 

where Aq is the amplitude, 6 and are the polar and azimuthal angles of the incoming light 
beam, and e^, are unit vectors along the Cartesian axis. For the setup of Fig. 5,^ = 0, and 
Eq. ^ reduces to 
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/ = e^/mc I Ao (cos 06^; + sin^ej^)- < ipjUpx + Py)|^i > P (19) 

i 

It turns out that the antibonding state a as well as the bonding state b is even and, therefore, 
recalling the aforementioned symmetry arguments, only the Px term contributes and Eq. ^ 
yields 

/ = eh/mcY, Bii\Ao < 4^f\px\ipi > P cos^ (f) (20) 

i 

which explains the shape of polar plots of Fig. 5 for states a and b; states c and c' must also 
be even in view of the shapes of the related polar plots. The preceding analysis is generalized 
straightforwardly to consider vector potentials with other non-zero components. If a non- 
zero Pz component is included, Eq. |20| is modified with a radially symmetric contribution. In 
this vein, if the mirror symmetry plane is rotated by 45° as is the case in YBCO (neglecting 
the effect of Cu-0 chains), that leads to polar plots which are rotated similarly. [p!8| , ^ 



Fig. 5 also shows the energy distribution curve (EDC) for 22 eV photons when ky is 
held fixed at M. The EDC for the case when the initial state width is taken to be a small 
constant is given by the thin line for reference in order to reveal the spectral peaks clearly. 
More realistically, the initial state width increases with increasing binding energy. This 
effect is modelled by the EDC of the thick line which shows that much of the structure for 
binding energies greater than ~ 100 meV is washed out by lifetime effects. 



D. Final state band structure, momentum matrix element 

Fig. 6 compares the one-step ARPES intensity at the M point for the antibonding state 
a with the corresponding result based on the bulk momentum matrix element of Eq. ^ 
and helps address a number of related issues. The small final state width used in these 
simulations allows the intrinsic spectral structure to be seen clearly. The final state bands 
shown as a function of k±_ in the lower portion of the figure provide a map of all available 
transitions since k± is not conserved in the photoemission process. [Note that the horizontal 
scale has been converted to refer to the photon energy.] It is striking that although these 
bands appear to form a semi-continuum, there are gaps in the underlying final state spectrum 
and the situation thus is quite unlike the free electron case [0. Also, k± dispersion is not 
always negligible with some bands dispersing as much as 0.5 eV despite the small size of the 
BZ in the z-direction. 

As already noted, it tends to be difficult to associate peaks in the one-step spectra 
with transitions to specific final states. The problem becomes exacerbated when states are 
broadened to incorporate finite lifetime effects. Considerable insight into the behavior of 
the ARPES matrix element can nevertheless be gained via the momentum matrix element 
< ijjf\p\ipi > of Eq. 0for bulk transitions. We emphasize however that < ipflpli'i > does 
not properly account for the processes of transport and ejection of the photoelectron as is the 
case in the one-step ARPES computation. Moreover, since < ipflplipi > involves the bulk 
wavefuctions ipi and ijjf, modifications of these wavefunctions in the vicinity of the surface 
and the effect of possible formation of surface states are not incorporated. Even so. Fig. 6 
shows that the one-step spectrum for the antibonding state a is rather similar to that based 
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on the momentum matrix element. Some differences in the intensities of certain peaks in the 
figure may be understood as consequences of the fact that the momentum matrix element 
has been evaluated at fc^ = while in the one-step result k± value is in general non-zero and 
varies from k± = to the zone boundaries at k± = ±7r/c for various peaks. Results for the 
bonding state b (not shown for brevity) are similar in that the level of agreement between 
the one-step and the momentum matrix element based spectra is comparable to that of Fig. 
6. 



V. NATURE OF TRANSITION MATRIX ELEMENTS AND ORIGIN OF 
SPECTRAL FEATURES IN TERMS OF CHARACTER OF INITIAL AND FINAL 

STATES 

A. Contributions of different atomic sites and related interferences 

Contributions M'^'^ of different atomic sites in the unit cell can be delineated via the 



decomposition of Eqs. |TJ-|16[ Table 1 and Fig. 7 provide illustrative examples of such an 
analysis. Considering x-polarized light, values of M]'^ for each of the 15 basis atoms for one 
particular transition from the states a and 6 at M are listed in Table 1. Fig. 7 on the other 
hand gives various significant contributions to all transitions from state a over the energy 
range of 10-27 eV. 

The strong influence of the mirror symmetry around the Ca layers is evident in Table 
1. For the particular transition considered, the magnitudes of MJ'^ are the same for the 
two mirror partners corresponding to the two 6 values, but the phases differ by vr. In fact, 
contributions cancel pairwise for the bonding state b, but interfere constructively for the 
antibonding state a. For many other transitions from a and b the situation is quite the 
opposite in that this sort of interference is destructive for a, but constructive for b. As we 
move away from M, these cancellation effects become less perfect as the phases do not always 
differ exactly by vr and/or the magnitudes of the two mirror partners are no longer the same. 

Figure 7 provides a somewhat broader perspective. Total matrix element \M^.\ obtained 
by summing over 7 and 6 is depicted pictorially for many transitions from a together with 
contributions from pairs of Cu, Ocu,x, Osi and Osr atoms; MJ for Ca, Sr, Bi and Ocu,y 
are less than 0.01 in magnitude and are not shown. Of a total of 64 possible transitions 
from state a over the energy range of Fig. 7, only 20 transitions that are shown possess a 
significant intensity. A few of the transitions are identified by numbers 1-7 for convenient 
reference. As already noted above in connection with Table 1, transition 1 is quite intense 
and arises almost completely from Ocu,x- Transition 4 contains an admixture of Ocu,x and 
Cu, while transition 3 contains additionally a contribution from Osi- Transition 5 contains 
contribution from all these atoms as well as Osr-, but as indicated by downarrows, Cu 
and Obi give negative contributions yielding a net reduction in the total intensity. These 
examples would make it clear that constructive or destructive interference can take place 
not only for a given atomic pair, but also between different atomic sites. We emphasize that 
the intensity variations in Fig. 7 reflect changes in the character of the final states since the 
initial state for all transitions is fixed to be the same. 
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B. Contributions of various angular momentum channels and interference effects 

therein 



Further insight can be gained by considering the angular momentum dependence of the 
matrix element, i.e. the quantities Mj*^ ^ of Eq. For illustrative purposes, we consider 
first the aforementioned transition 1 with the help of Tables II-V, followed by a few comments 
concerning the behavior of transitions more generally via Fig. 8. 

The weights in different angular momentum channels (defined as \C]^\'^) for the initial 

5 

states a and h associated with various atomic sites are given in Tables II and III, while the 
final state for transition 1 is characterized similarly in Table IV EBI]. We emphasize that 



these weights depend upon details of how space within the unit cell is ascribed to different 
atoms and the size and shape of the interstitial region. For this reason, we should keep in 
mind that the results of Tables II-IV possess inherent uncertainty, even though these are 
representative of the character of the states involved. As expected, wavefunctions of both 
a and h are seen from Tables II and III to be predominantly Cu{dx'i-y2) and Ocu,x{Px) like; 
Obi as well as Osr states hybridize more strongly with the bonding level h compared to a. 
The final state in Table IV continues to be weighted in favor of Cu and Ocu,xi but being 17.7 
eV above the initial states, there is substantial change with respect to distribution among 
angular momentum channels and atomic sites. 

We turn next to Table V and the matrix element for transition 1 from state a. Of the 
large number of possible terms (~ 300) under the summations in Eq. only the few 
elements Mj^^^, listed in Table V possess substantial magnitude. The remaining elements 
are either zero due to selection rules or are negligibly small. The px dx'^-y^ channel 
on Ocu,x is seen to be the strongest with a contribution (in arbitrary units) of 0.27. It is 
most striking that the net contribution from Cu is only 0.02 even though the initial state is 
dominated by Cu. This is the result of strong interference effects: On Cu, different angular 
momentum channels interfere destructively as manifest in the positive and negative signs of 
elements, while on Ocu,x, there is a constructive interference yielding a large total intensity 
of 0.43. These observations hint that the ARPES matrix element can display remarkable 
site selectivity properties. For transition 1 under discussion, for example, almost all the 
photoelectrons emanate from Ocu,x sites with little contribution from other sites in the 
crystal due to complicated interferences. 

Fig. 8 considers all transitions from a over the energy range of 10-27 eV. The total 
matrix element \Mx\ and the portion \M]\ associated with Ocu,x in the top two panels of 
the figure are reproduced from Fig. 7 for ease of discussion; the total Ocu,x intensity is 
decomposed into angular momentum channels in the three lower panels. By comparing the 
lengths of arrows in Figs. 8(a) and 8(b), we see that for most transitions (but not all) Ocu,x 
component is strong. The effect however is not as dominant as in the case of transition 1; 
for example, only about 60 percent of transitions 2-4 come from Ocu,x and transition 6 has 
essentially a zero Ocu,x part. Moreover, even though different angular momentum channels 
interfere constructively for most transitions in that arrows are mostly pointed up in Figs. 
8(c)-(e), this is not always the case as some arrows are pointed downwards. In particular, 
for transition 5 roughly half of the intensity in s channel is destroyed by inteference 

from px dx2-y2 channel, and for transition 7, these two channels destroy each other nearly 
completely. We note finally that, of the many factors involved in defining M'^^^, on the 
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right side of Eq. 17, much of the energy dependence depicted in Fig. 8 originates from 
the behavior of the wavefunction coefficients of the final states (C^?); the Gaunt coefficients 
Gl,u '^^i^c constants and the quantities B'^^^ involving integrals of radial wavefunctions are 
slowly varying functions of energy. 

C. k± dependence 

So far we have focused on the excitation of states a and h at the M(7r,0,0) point, i.e. 
for A;_|_ = 0. Changes in the momentum matrix element as a function of kx_ for several 
selected transitions from the initial state a are considered in Fig. 9. Transitions 1-4 have 
been discussed previously in connection with Figs. 7 and 8. Transitions 2b and 3b however, 
which lie quite close in energy to transitions 2 and 3, respectively, were not shown before 
since their intensity is zero or near-zero at = 0, but as Fig. 9 shows, this is not the case 
for 7^ 0. Transitions 1 and 4 decrease slightly in intensity as fc^ increases from zero to 
tt/c. On the other hand, transitions 2, 2b, 3 or 3b undergo large changes. Observable effects 
after spectra are smoothed with realistic broadenings of the order of a few eV's will likely 
be much smaller since the intensities of transitions 2 and 2b (or similarly 3 and 3b) taken 
together are expected to change little with k^. 

We have carried out additional computations at k-points on the F — M line away from 
M . The results continue to display complicated interference effects involving pairs of atoms 
and/or various atomic sites, essentially along the lines of the preceding commentary con- 
cerning emissions from M, and require little further elaboration. 

VI. SUMMARY AND CONCLUSIONS 

We have investigated the nature of the ARPES spectrum from BISCO using the first- 
principles band theory framework. Our focus is on understanding emission from electronic 
states along the F — M symmetry line and particularly from the bonding and antibonding 
combinations of bands formed from the states associated with the Cu02 bilayers which have 
been the subject of considerable recent discussion and controversy. Salient features of the 
ARPES spectrum delineated are as follows. 

1. The intensity of the antibonding as well as the bonding state at M displays a strong 
dependence on photon energy even though the two states differ significantly from each other 
in this respect. This result indicates that the bonding state will be highlighted at certain 
photon energies while the antibonding state will be intense at other energies; moreover, at 
photon energies where emission from either state is relatively weak, we would expect the 
effects of deviations of the underlying lattice from perfect tetragonal symmetry (i.e. via 
various modulations) to become more visible. 

2. The spectra depend strongly on the polarization of the incident light. At the M 
point, for states in the vicinity of the Ep, the intensity is maximum when light is polarized 
along the x- or the Cu — bond direction, but becomes nearly zero upon a 90° rotation 
of the polarization vector perpendicular to the plane of the detector. Simple arguments are 
presented to explain the cos^0 (where is the angle between the x-axis and the polarization 
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vector) variation observed in the one-step computations which yield a polar intensity plot 
in the shape of the number eight in BISCO and other cupratcs. 

3. An examination of the final states indicates that, despite the appearance of a semi- 
continuum, there are gaps in the final state spectrum and that the k± dispersion of some 
bands is as large as 0.5 eV even though the Brillouin zone size is small along the z-direction. 
The situation thus is quite unlike the free electron case. 

4. We show how considerable insight into the nature of the specific ARPES excitations 
can be adduced in terms of the momentum matrix element for bulk transitions within the 
solid in combination with constraints appropriate for a particular ARPES measurement. 
Such transition specific information tends to be difficult to ascertain within the more com- 
prehensive one-step framework, especially in the presence of finite initial and final state 
dampings. 

5. With the preceding point in mind, we have studied the nature of the momentum 
matrix element for emission from the CuOi plane related antibonding state a and the bond- 
ing state 6 at M at length over the transition energy range of 10-27 eV. We show how the 
matrix element, |Mp = | < ■?/'j|p|-?/'/ > p, can be decomposed usefully into contributions 
from various atomic "pairs" in the unit cell (i.e. an atom and its mirror partner obtained 
by reflection through the Ca layer about which the BISCO lattice is symmetric), and fur- 
ther into different angular momentum channels. Contributions from the two atoms in the 
aforementioned atomic pairs display remarkable effects of interfering constructively for some 
transitions and destructively for others. The net contribution from the Ocu,x pairs seems 
to dominate the spectrum in many cases even though the initial and final states involved 
possess large weights on the Cu sites. These strong interference effects continue to be present 
in different angular momentum channels. The analysis indicates that many transitions are 
controlled essentially by just the single — > d^2_y2 channel on the Ocu,x sites. This striking 
result hints that the ARPES matrix element may possess most interesting site selectivity 
properties. 

6. The intensity of many transitions increases (or decreases) siibstantially as a function 
of k±. However, these variations are usually coupled with nearby transitions whose intensity 
decreases (or increases) so that the net contribution is generally not expected to show a 
strong k± dependence. 

In summary, the present analysis makes it clear that the complicated nature of the 
ARPES matrix element in complex materials hides considerable potential for extracting 
useful information via this spectroscopy about the nature of specific states. Studies of other 
cuprates along these lines should allow an identification of robust aspects of the electronic 
structure related to the Cu02 planes. 
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TABLES 



TABLE L Magnitude and phase of the matrix element MJ' is given for various sites 7 and 6 
for transition 1 in Fig. 7 (or Fig. 8) for both the antibonding initial state a and the bonding initial 
state b at M. Note that Figs. 7 and 8 only refer to state a, but that this table also gives results for 
the state b. 
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TABLE II. Wavefunction of the antibonding state a at M-point is characterized in terms of 

weights (in percent) defined as J^lCjfl"^ associated with different sites 7 (see Eq. |l^. Not all 

s 

angular momentum channels are listed. Channels with small weights are omitted for brevity; the 
total weight for this reason adds to less than 100 percent. [46| 



7 


Tot. 


s 


Pz 


Px 




dx^—y'^ 


Ca 














Bi 


12 




7 


5 






Sr 














Cu 


43 


3 






12 


28 


Ocu,x 


24 






24 






Oc'u,y 














0B^ 


4 




4 








Osr 


6 




6 









18 



TABLE III. Same as the caption to Table II, except that this table refers to the bonding state 
b at the M-point. 
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TABLE IV. Same as the caption to Table II, except that this table refers to the final state of 
transition 1 in Fig. 7 (or Fig. 8) involving the excitation of the antibonding initial state a. 
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TABLE V. Components J2^2ll' momentum matrix element (see Eqs. 14-17) for 

s ' ' 

transition 1 marked in Fig. 7 (or Fig. 8) which involves excitation of the antibonding initial state 
a at M-point. Only sites and channels with significant weight are listed. 
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FIGURES 



FIG. 1. Body-centered tetragonal crystal structure of Bi2212. (a) and (b) give (110) and (100) 
cross-sections through the structure, respectively, (c) depicts the arrangement of atoms in different 
layers. Various basis atoms are identified by the pairs of indices (7, S) discussed in the text. 



FIG. 2. Band structure of Bi2212 in the vicinity of the Fermi energy Ep (shown by the 
horizontal dashed line) along the symmetry lines T — M — X — T. Note that the Bi-0 pockets 
around M in the standard LDA band structure have been lifted above the Ep; see text for details 
of crystal potential, a and b are the antibonding and bonding states at M, respectively, arising 
from the Cu02 planes, c and c' are Bi-0 related states. 



FIG. 3. ARPES intensity for emission from the Cu02 plane related antibonding and bonding 
states at M function of the photon energy. 



FIG. 4. ARPES intensity from the antibonding state a at M as a function of the photon 
energy for four different values of the width of the final state determined by the parameter |41 



FIG. 5. Polarization dependence of the ARPES spectra at M for 22 eV light. Energy distri- 
bution curve (EDO) including realistic initial state broadening for light polarized along the x-axis 
is shown (thick line); thin line gives the EDC computed by using an artificially small damping of 
initial states to reveal the spectral peaks a, 6, c and c' clearly. EDO's have not been convoluted 
with the Fermi function. Inserts in the shapes of number eight on top of various peaks are polar 
plots of peak intensity as the direction of the polarization vector (1^) is varied and possess the cos^cf) 



dependence in accord with Eq. 2C. </> is measured with respect to the x-axis or the Cu — Ox bond 
direction; as seen in the inset in the top part of the figure, the detector (D) is assumed to lie in 
the vertical plane through the x-axis. 



FIG. 6. Computed one-step ARPES intensity in BISCO at M from the antibonding state a as 
a function of the photon energy is compared with the corresponding result based on the momentum 
matrix element |Mi,p for bulk transitions (see Eq. 14). A small final state damping is used to 



reveal features clearly. The relevant final state bands are shown as a function of k_\_ and energy in 
the lower figure where the energy scale has been shifted to correspond to photon energy. 
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FIG. 7. Total matrix element \Mx\ (topmost panel) for various transitions from the antibond- 
ing state a over the transition energy range of 10-27 eV is decomposed into contributions from 
various "pairs" of atomic sites 7 (panels (b)-(e)) as described in Eqs. 14-^. Down arrows indicate 
a negative value of MJ. Sites not shown yield a negligible contribution. Transitions labelled 1-7 
(connected by dashed lines for ease of identification across the panels) are discussed in the text to 
highlight various interference effects. 



FIG. 8. Same as the caption to Figure 7 with the two top panels (a) and (b) reproduced from 
Figure 7 for ease of reference. Panels (c)-(e) give the decomposition of Ocu,x results in (b) into 
various angular momentum channels J2 -^J l L' (^^^ ^1- HI) • 



FIG. 9. k± dependence of the total momentum matrix element \Mx\ for a few selected transi- 
tions from the antibonding initial state a at M. Transitions 1-4 are the same as those in Figs. 7 
and 8. Transition 2b (or 3b) possesses energy close to 2 (or 3) but its intensity is zero (or quite 
small) for k± = 0, or equivalently, at M(7r,0,0). 
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